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Abstract 

The space-time evolution of the (l+l)-dimensional viscous hydrodynamics with an initial quark- 
gluon plasma (QGP) produced in ultrarelativistic heavy ion collisions is studied numerically. The 
particle-emitting sources undergo a crossover transition from the QGP to hadronic gas. We take 
into account a usual shear viscosity for the strongly coupled QGP as well as the bulk viscosity 
which increases significantly in the crossover region. The two-pion Hanbury-Brown-Twiss (HBT) 
interferometry for the viscous hydrodynamic sources is performed. The HBT analyses indicate 
that the viscosity effect on the two-pion HBT results is small if only the shear viscosity is taken 
into consideration in the calculations. The bulk viscosity leads to a larger transverse freeze-out 
configuration of the pion-emitting sources, and thus increases the transverse HBT radii. The results 
of the longitudinal HBT radius for the source with Bjorken longitudinal scaling are consistent with 
the experimental data. 
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I. INTRODUCTION 



The experimental results of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory (BNL) indicate that the matter produced in the central collisions of 
Au + Au at ^snn =130 and 200 GeV is a strongly coupled quark-gluon plasma (sQGP), 
and behaves like a perfect liquid Recently, the studies of the dissipative fluid dynamics 
become very hot in high energy heavy ion collisions 2|-|20j. because people want to know 
why the sQGP exhibits an almost perfect fluid property and how to probe the viscosity 
effects on experimental observables. 

The relativistic formulism of dissipative fluids was originally derived by Eckart 2l| and 



Landau and Lifshitz 



22] . Their theories contain only the first-order terms of the dissipative 



quantities and therefore are referred to as the first-order theories of dissipative fluids. The 
theories of dissipative fluids which include the terms up to the second-order of the dissi- 
pative quantities were developed by Miiller (non-relativistic) [23[ and Israel and Stewart 
(relativistic) 24j. These second-order theories can avoid the problem that the first-order 
theories may not be satisfied causality sometimes. So they are also called the causal theories 
of dissipative fluids dynamics. 

Although the Israel- Stewart second-order formalism for relativistic dissipative fluids 24 1 
was established thirty years ago (in the 1970s), its numerical implementation in high energy 
heavy ion collisions begins in the 21st century. In Ref. {25 1, the viscous hydrodynamics 
based on the second-order theory is first used to study the system expansion in relativistic 



heavy ion collisions with the Bjorken scaling hypothesis 



261 ] . Recently, the causal dissipative 



fluid dynamics hasjpeen used in high heavy ion collisions for the investigations of the effects 

4l6| and both shear and bulk viscositie s |17H 20| on the transverse 
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of shear viscosity 
momentum distribution |4H7 
Twiss (HBT) interferometry 



Jll 16tj2l 



J20J], elliptic flow 



19], and Hanbury-Brown- 



As is well known, a dissipative fluid has not only shear viscosity but also bulk viscosity. 
Recent researches jl?! \v\ indicate that the effects of the shear and bulk viscosities on the 
elliptic flow are completely different in the RHIC heavy ion collisions. The shear visc osity 
suppresses the V2 at larger transverse momentum but the bulk viscosity increases it 17J, Il8 |. 
Unlike elliptic flow, which reflects the anisotropic pressure of the system, HBT correlations 
are related to the space-time structure of the particle-emitting source. The investigations 



2 



of the shear and bulk viscosity effects on HBT interferometry may probably provide the 
messages of the shear and bulk viscosities from another aspect. 

In the present work, we begin with the equations based on the Israel-Stewart theory for 
the net baryon free system formed in ultrarelativistic heavy ion collisions. For the cen- 
tral collisions, it is assumed that the system has the Bjorken cylinder geometry (expanding 



isotropically in transverse and satisfying longitudinal boost-invariance) [27H29I] and under- 
goes the crossover transition from the QGP to hadronic gas. We shall take into account 
both the shear and bulk viscosities and solve the viscous hydrodynamic equations by the 
relativistic Harten-Lax-Leer-Einfeldt (RHLLE) algorithm 29|-|32j. By applying a simulated 



analysis of two-pion HBT interferometry [32l . 



33] to the (l+l)-dimensional viscous hydro- 



dynamic sources, we investi gat e the effects of the shear and bulk viscosities on the HBT 



radii R on t, Rside, and -Ri ong 
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351 ] . The HBT analyses indicate that the viscosity effect 



on the HBT results is small if one considers only the shear viscosity. The bulk viscosity 
leads to a larger transverse freeze-out configuration of the pion sources, and thus increases 
the transverse HBT radii. The results of the longitudinal HBT radius for the source with 
Bjorken longitudinal scaling are consistent with the experimental data. 

This paper is organized as follows. In Sec. Ull we review briefly the evolution equations 
of the relativistic dissipative fluid dynamics for the net baryon free system with the Bjorken 
cylinder geometry. We present the expressions of the relaxation equations for the shear 
tensor and bulk pressure for the system considered. The detailed derivations are given 
in Appendix A. In Sec. 1111} we discuss the numerical algorithm for solving the viscous 
hydrodynamics. The equations of state and initial conditions used in our calculations are 
outlined. We present our numerical results of the viscous and ideal fluid evolutions. In Sec. 
\IV\ we perform the two-pion HBT interferometry for the viscous hydrodynamic sources. We 
investigate the viscosity effects on the HBT radii for the viscous hydrodynamic sources with 
only shear viscosity and with both shear and bulk viscosities. We also investigate the HBT 
results for the viscous sources with lower and higher freeze-out temperatures and initial 
energy densities. Finally, a summary and discussions are given in Sec. [V] 
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II. EQUATIONS OF RELATIVISTIC DISSIPATIVE FLUID DYNAMICS 



First, we briefly review the evolution equations of relativistic dissipative fluid dynamics 
for the system with the Bj or ken-cylinder geometry and zero net baryon density The basic 



equations can be referred to 



-4 



36). 



The evolution equations of hydrodynamics come from physics conservations. For the 
baryon free system, we consider only the energy-momentum conservation in the hydro- 
dynamics. It is convenient to describe the Bjorken- cylinder system with the coordinates 
(r, p, (ft, 77), where r = \/t 2 —z 2 is the longitudinal proper time, p and <fi are the polar coor- 
dinates in the plane transverse to the beam direction z, and r\ = | ln[(t+z)/(t— z)] is the 
space-time rapidity. The metric tensors g^ v and are expressed in the frame of these 
curvilinear space-time coordinates as 



gT = diag(l,-l,-l/p 2 ,-l/r 2 ) : 
= diag (l, -1, -p 2 , -r 2 ) . 



(1) 



Some notations used in this paper are 



r 7 

1 a/3 



^9 la {d a g/3a + dpg^ - d a g ai3 ). 



= dX» 



D = u>"d 



A 4 * 



ee - vPu", 

\7 <tA u u> ee V { ^u u) - -A^d x u } 
3 

A ( a p) _ + A f>«y 2i 



(2) 

(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



were denotes the covariant derivative, d a = d/dx a and = 7(1, v) are the 4-derivative 
and 4- velocity (7 = l/y/i — v 2 ), is the Christoffel symbol, and A^ a ^ denotes the sym- 
metrization of the quantity A a/3 . 



We adopt the Landau and Lifshitz frame 



221 ] in which the local energy flow is zero. The 



energy-momentum tensor of a fluid cell in relativistic dissipative hydrodynamics is 



T"" = evPu" - (p + II) A"" + tt^, 



(9) 
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(11) 



where, e is the energy density, p is the local isotropic pressure, IT is the bulk viscosity 
pressure, and is the shear stress tensor. Using the covariant derivative form, the energy- 
momentum conservation can be reexpressed as 

dpTi™ = dpT 1 *" + T^ x T Xu + T^T^ = 0. (10) 

For the Bjorken cylinder, the longitudinal velocity is given by v z = z/t and the hydrody- 
namical solution at an arbitrary longitudinal coordinate z can be obtained by the Lorentz 
boost with the rapidity rj = tanh _1 (z/t) from the system's transverse evolution at z = 
27H29I. 1371] . So we need only to solve the hydrodynamical equations at z = (rj = 0) in this 
case. Considering u v = at z = and assuming i$ = 0, Eq. ffTUj) becomes 

d T T TT + d p TP T + T T m T^ + T^ T T TT + rJ p T" r = 0, 
d T T Tp + d p TPP + T P H T^ + Y\ p TPP + r^T r " = 0. 
Introducing the quantities 

p p = p + n + 7r w / 7 2 , (12) 

^=p + II + 7r*Y, (13) 
P^p + n + vr^r 2 , (14) 

E = (e + P p )7 2 -P P , (15) 
M=(E + P> p , (16) 
the non zero components of the energy-momentum tensor can be expressed as 

T rr = E ^ T rp = rppr = ^ ^ 

T p P = Mv p + V p , (18) 

j,(j>(j> __ rprm _ !_± Mg\ 

and, Eq. ljTT]) can be rewritten as 



d T £ + d p [(E + V p ) v p ] = —(E + V p ) + -) + Up (20) 

V p T / T 

d T M + d p (Mv p + V p ) = -m(^ + -)- ( V p- V *\ (21) 

V p T J p 

The coupled equations ( 120!) and (I2T]) are the evolution equations of the dissipative hydro- 
dynamics for the Bjorken cylinder system, which can be numerically solved by the RHLLE 



algorithms [29l-l32j . The quantities V p , V^, and V n include the effects of dissipation and they 
reduce to the isotropic pressure p for perfect fluids. 

Next, we present the relaxation equations that the dissipation quantities (ir^ and II) are 
satisfied for our system. The detailed derivations of the relaxation equations are given in 
Appendix A. In order to obtain the system evolution we need to solve the equations ( 120]) 
and (12 ip together with the relaxation equations. 

In the Landau and Lifshitz frame and for the baryon free system, the dissipative effect 



from the heat conduction can be neglected [2 
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121 ] . On the other hand, one can also 



ignore the effect of vorticity because it is small for the longitudinally boost-invariant system 
. . n this case, the relaxation equations for the shear tensor tc pu and bulk pressure 
n can be written as Qfl S, 36] 

^i"'T>T,lA — 



2 V T) 1 



(22) 



1„~„, r r n u } 



mDTl + n = -C dX ~ 2 n C Td x { 



(23) 



where , rj and ( denote the shear and bulk viscous coefficients, r n = 2f3 2 rj and r n = (3 ( 
are the corresponding relaxation times, 02 and /3q are the relaxation coefficients for the 
expansion of entropy flow up to the second-order of the dissipative quantities, and T is the 
temperature. 

For the symmetric Bjorken cylinder and at z = (vP = u v = 0), Eqs. (122]) and (1231) can 
be written in the forms which are suitable for solving numerically as (for detailed derivations, 
see Appendix A) 



d 00 8t pp 1 

7T r + v p^t~ = — 1 

or op "yr n 



-T PP + 2?)(J PP ) - 



T 



PP 



27 



e 



rjT 



(24) 



At* + v ^ 
dr p dp 



l-T*+ + 2rio**) 

ir n ' 27 



T n \7)T 



(25) 



d dr m 
_ T nn _|_ v 

dr 9 dp 



1fT„ 



'-T vr, + 2?ja 



T 



1 in 



2 7 



rjT 



(26) 



d_ 
dr' 



on 



— n+ vF = — (-n 

dp 7T n 



CO) - 1 

2 7 



rn V CT 



(27) 
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where r pp , r^, and r m are three introduced quantities whith the relations to the nonzero 
components of 7T flu as 

\ 



/ ~f 2 v 2 p r pp -i 2 v p tPP 



71 



7 2 V p T PP 



















2^r/ 



T T 



(28) 



V 



where 



T' 



PP 



'/'/ 



from the traceless condition ir^ 

p 







e 



d(7 ^p) 
<9r <9p 



25|]. In Eqs. ([1 



PP 



-e 



p 



i 

—i 

3 



a 



vv 



T 



7 1 

l + -e 



(29) 

(30) 

(31) 
(32) 

(33) 



and D = u^d^ reduce to 7(<9 r + v p d p ). 

When performing numerical calculations we need to solve the evolution equations (|20|) 



and (Jzlj) and the relaxation equations (J24T) - (j27j) simultaneous 
of the shear viscosity to the entropy density satisfies rj/s>l/A7i 



38 



or the sQGP, the ratio 



401 ]. and the ratio of the 



bulk to shear viscosities is proportional to the deviation of the sound velocity square from 



18 



41 



4311 . Based on the calculations of 



that of the ideal hadronic gas, C/v — ~ K ( v s ~ 1/3) 
the M = 2* gauge theory, the value of k is between 3.142 and 4.935 43]. In our calculations 
we take rj/s = 2C s (l/47r) (C s = 1, 2) and k 

bulk viscosities are taken and Tn = r n as in Refs. (7j,Hl|, 

simplicity. 



4.75. The relaxation times for the shear and 

nn 
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16 



18 



19 



24| for 



III. NUMERICAL SOLUTION OF VISCOUS HYDRODYNAMICS 

When solving the viscous hydrodynamic equations numerically, we encounter the follow- 
ing two types of equations: 



d T U + d p F{U) = G{U) 



(34) 



and 



d T U + v p d p F(U) = G(U). 



(35) 



For the evolution equations (!20|) and (12TT) . U in Eq. (134")) can be replaced by E or M, and for 



the relaxation equations (jUj) - (1271) . U in Eq. (135]) can be rep_ 



aced by r w , r^, r w , or n. Eq. 



(134p can be solved directly by the RHLLE algorithm 29h32| . For Eq. (1351) we can change 
it to the type of Eq. ( |34"1) by adding the term F(U) d p v p to both sides of the equation, and 
then solve it with the RHLLE algorithm. In the calculations, we solve simultaneously the 
six coupled equations of (j20) . (j2Ij). and (12D - (f2TI) for the quantities (E, M, r rr , r^, r^, n) 
at each time step. The width of time step is taken to be At = 0.04 fm/c and the width of 
space step is taken to be 0.99 At fm. The set of the coupled equations (12"0|) . (12"TT) . and (12"4"|) 
- ( 1271) is closed by the equation of state (EOS), and the transverse velocity and local energy 
density satisfy 

MM 2 
v p = — — , e = E — — — . (36) 

p e + p; e + p p y 1 

The EOS used in the calculations is the parametric EOS which combines hadron resonance 



gas at low temperatures with lattice QCD at high temperatures |44|. Figured] (a) shows 
the energy density e and pressure p as functions of temperature. Figure [I] (b) shows the 
square of sound velocity v 2 and the difference between e and 3p, A = e — 3p, as functions of 
temperature. The transition temperature T c is taken to be 170 MeV. 



On the basis of the Bjorken picture 26], the system evolves hydrodynamically after a 
proper time To. In our calculations, we take the initial time To = 0.6 fm/c, and use a 
Gaussian initial transverse distribution of energy density at z — 0, 

e(p,z = 0)=e eM-p 2 /2a ), (37) 

where Eq is taken to be 30 GeV/fm 3 and <Jq is taken to be 3.3 fm. The initial shear tensor 
components are taken to be r pp (p) = r^(p) = and r w (p) = as in Ref. ^|. We 

find that the system evolution is almost independent of the initial value of the bulk pressure 
because it is almost zero at the initial temperature much higher than T c . So we take the 
initial bulk pressure to be zero in the calculations. 

Figures [2] (a) and (b) show the temperature and transverse velocity profiles in z = 
plane for the viscous fluid with only shear viscosity (C s = 2). Correspondingly, Figs. [2] (c) 
and (d) show the temperature and transverse velocity profiles in z = plane for the fluid 



S 




FIG. 1: (Color online) (a) The energy density e and pressure p as functions of temperature, (b) 
The square of sound velocity v$ and the difference A = e — 3p as functions of temperature. 
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FIG. 2: (Color online) The temperatures and transverse velocities of the viscous fluids with only 
shear viscosity (C s = 2) and with the shear and bulk viscosities, in z = plane. The labeled time 
is respectively t = r , t n = r + 2 n (20Ar) (n = 1, 2, 3, 4, At = 0.04 fm/c). 
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FIG. 3: (Color online) The ratios of t pp , t^, t vv , and II to £o as functions of the transverse 
coordinate and time for the viscous fluid (C s = 2). Here, the time is taken to be tq (solid black), 
r + 6 At (dashed black), r + 12Ar (solid red), r + 24Ar (dashed red), r + 48Ar (solid blue), 
and ro + 96Ar (dashed blue), respectively. (Ar = 0.04 fm/c) 

with the shear and bulk viscosities. In Fig. [2j the dashed lines are for the ideal fluid for 
comparison, and the time to = T o, t n = r + 2"(20Ar) (n = 1,2,3,4, Ar = 0.04 fm/c). 
For the only shear viscosity case the viscous fluids cool a little more slowly than the ideal 
fluid for small time, and the transverse velocities of the viscous fluid are higher than the 
corresponding velocities of ideal fluid. Once the bulk viscosity is taken into account, one 
can see that the viscous fluid cools more slowly and its transverse velocities are lower than 
the corresponding velocities of ideal fluid. This is because that the bulk viscosity decreases 
the system pressure (see Eq. (EJ) and notice II < 0) and therefore decreases the gradient of 
the pressures in and out of the system. 

In Fig. [3], we show the ratios of the viscosity quantities r pp , r^, t vv , and II to the central 
initial energy density Eq, as functions of the transverse coordinate and time for the viscous 
fluid (C 8 = 2). In FigJU we draw the isotherms for the viscous (C s = 2) and ideal fluids in 
z = plane. It can be seen that the contours for the only shear viscous fluid are slightly 
larger than that of the ideal fluid, and the bulk viscosity lets the contours of the viscous 
fluid be much larger than those of the ideal fluid. This is consistent with the results in Fig. 



After knowing the transverse evolution at z — (t — r), we can obtain the temperature 
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FIG. 4: (Color online) The isotherms for 0.7 and 0.8 T c in z = plane. 



and velocity at arbitrary z for the Bjorken cylinder system by 27H29I. l37l | 

T{t,p,z) = T{r,p,0), 



T 



v p (t,p,z) = v p (r,p,0)-, v z =-. 
Then, we can calculate the two-pion HBT correlation functions by a simulation 



(38) 
(39) 
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IV. PION INTERFEROMETRY ANALYSIS 

The two-pion HBT correlation function is defined as the ratio of the two-pion momen- 
tum distribution P{k\, k 2 ) to the the product of the single-pion momentum distribution 
P(fci)P(/c 2 ). For a chaotic pion-emitting source, P(ki) (i = 1,2), and P(ki,k 2 ) can be 
expressed as [45] 



p(h) = Y,A\h,x t 



(40) 



x, 



P(h,k 2 )= f(h,k 2 ;X 1 ,X 2 ) 

Xl,X2 



(41) 



where A(ki,Xi) is the magnitude of the amplitude for emitting a pion with 4-momentum 
ki = (ki,Ei) at Xi and $(&i, k 2 ; Xi, X 2 ) is the two-pion wave function. Assuming that the 



11 



emitted pions propagate as free particles, k 2 ; Xi, X 2 ) is simply 



$(k 1 ,k 2 ;X 1 ,X 2 
1 

+ A(ki, X 2 )A{k 2 , Xl )e lkl - X2+ik2 - Xl 



A(ki, Xx)A(k 2 , x 2 )e lkl - Xl+lk2 - X2 



(42) 



For a set of variables of the relative momentum, {qj} (qj = \hi — h 2 \j) , the two-pion 
correlation function can be expressed as 33] 

Cor({ qj }) 



c({ qj }) 



where 



and 



Uncor({gj}) ' 

Cor({g,})= f dk x dk 2 P{k h k 2 )\[K\^-Hj-^) 

J 3 

Uncor({ ? ,}) = J dkidk a P(A;i)P(fe) n<J(|k! - k 2 |, - 



(43) 



(44) 



(45) 



are the correlated and uncorrected pion pair distributions with {qj}. 

In our simulated calculations, we first generate the pion momentum k on the freeze-out 
hypersurface S(X) with temperature Tf, according to the probability of the Cooper- Frye 
formula [46| 

P(k)<x J (46) 

For the viscous fluid, the distribution f = fo + 5f. We take fo for the ideal fluid as the 
Boltzmann distribution, and take Sf for the shear and bulk viscosities as 10] 



(47) 



7T 



fii/ 



7] — t' 1 ' 1 T ppr ) 2 v p cos cosh i] r pp/ y 2 v p sin cosh 7] IP* 7 cosh 77 sinh 77 ' 



In Cartesian frame, Tt pv can be expressed in terms of the cylindrical variables as 

/ IF^cosh 2 " - '~/ v ! ~ 

T pp/ y 2 v p cos cosh 77 IP* cos 2 + r** IP* cos sin r ppr y 2 v p cos sinh 77 
r pp 7 2 f p sin cosh 77 IP* cos sin IP* sin 2 + r** T ppr y 2 v p sin sinh 77 
y n pr? cosh 77 sinh 77 r PPr y 2 v p cos sinh 77 r PP/ y 2 v p sin sinh 77 II P?? cosh 2 77 — r pp/ y 2 v 2 J 

(48) 

where fi,v = t, x, y, z, IP 7 ? = t pp j 2 v 2 + r w , and IP* = r w 7 2 - r**. 
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FIG. 5: (Color online) The two-pion HBT results for the viscous (C s = 2) and ideal hydrodynamic 



sources with Tf = 0.7T C and Tf = 0.8T C . The experimental data are from PHENIX 
481 ] Collaborations. 
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After obtaining the momenta fcj 



1,2) of the pions emitted on the freeze-out hy- 



per sur 
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ace, we can construct the correlation functions for the variables g out , q S ide, and qi ong 
35|, C(g out , (/side, Qiong), based on Eqs. fHOj) — f|45|) . by summing over ki and k 2 for each 

(^out) Qside; (Zlong 

correlation function with the parametrized formula 



bin. Then, we extract the HBT radii R Q ut, Rside, and i?i ong by fitting the 



C(q 

out; Qside; ^long 



l + \e- q ° 



2 p2 _„2 r>2 _„2 d2 
ut-n-out "We^side Vicing ^long 



(49) 



Figure |5] shows the HBT radii R Q ut, -Rside, -Riong, an d A parameter as functions of the 
transverse momentum of the pion pair, kx = |krx + ^-2t\/2, for the viscous and ideal hydro- 
dynamic sources with the freeze-out temperatures 0.7 and 0.8T C . For the viscous sources, the 
parameter C s is taken to be 2. For comparison the RHIC experimental HBT results (43, |3| 
are also shown in the figure. It can be seen that the HBT radii for the lower freeze-out 
temperature are larger than those for the higher Tf. The viscosity increases the transverse 
HBT radii R OVLt and -R S ide for the the source with the lower freeze-out temperature. By 
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comparing with the experimental results, we find that the results of the longitudinal HBT 
radius -Ri on g for the hydrodynamic sources with Tf = 0.8T C agree well with the experimental 
data. It indicates that the Bjorken scaling hypothesis is suitable in the description of the 
source longitudinal evolution. In our calculations, the width of the Bjorken rapidity plateau 
is taken to be 1.3. On can also see that the transverse HBT radii i? Q ut and -R S ide as functions 
of fcr are flat for both of the two freeze-out temperatures. However, the ratio -R ut/-Rsido 
decreases to about one for the higher Tf. We did not consider the source coherence and the 
Coulomb interaction between the pions, so the A results of the hydrodynamical sources are 
about unit. 

In Table I we list the HBT results for the ideal and viscous hydrodynamical sources with 
Tf = 0.8T C . Here the results of the 'vise fluid 1' and 'vise fluid 2' are for the viscous sources 
with only shear viscosity for C s — 1 and C s = 2. The results of the 'vise fluid 3' and 'vise 
fluid 4' are for the viscous sources with both the shear and bulk viscosities for C s = 1 and 
C s = 2. One can see that for the higher freeze-out temperature, the viscosity effect on the 
HBT results for the only shear viscous sources are negligible, and the viscosity effect for the 
source with both the shear and bulk viscosities may let the HBT radii be a little larger than 
those of the ideal fluid source. 

Further, we examine the HBT radii for the viscous hydrodynamical sources with lower ini- 
tial energy density (eq = 20 GeV/fm 3 ) and higher initial energy density (eq = 40 GeV/fm 3 ). 
The results are shown in Fig. El where the freeze-out temperature is 0.8T C . It can be seen 
that the results of the transverse HBT radius -R S id e for the higher initial energy density are 
larger than those for the lower initial energy density. This leads to the smaller values of 
the ratio -R ut/-Rside for the source with higher e . The longitudinal HBT radius -Ri on g also 
increases with the initial energy density. However, the transverse HBT radii as functions of 
the transverse momentum kx are more flat than those of the experimental data, both for 
the lower and higher initial energy density sources. 



V. SUMMERY AND DISCUSSIONS 



Based on the Israel-Stewart second-order theory of relativistic dissipative fluid dynamics 



241 ]. we investigate the space-time evolution of the viscous hydrodynamics with the Bjorken 



cylinder geometry. The EOS used in the hydrodynamic calculations combines hadron res- 
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TABLE I: The HBT results for the ideal and viscous hydrodynamic sources with Tf = 0.8T C . 



Kt {(je V /cj 


sources 




-Ksidc (nn) 


-Klong (tm) 


-"out / -n-side 


A 


U.Z4 


ideal fluid 


a ko _i_ n o/i 

4.D» ± U.U4 


r m J. n n/i 
O.Ul ± U.U4 


k on 4- n ok 
o.oy ± U.Ud 


n O/l 4- O OO 

u.y4 ± U.Uz 


n oi 4- n m 
u.yi ± U.U1 




vise fluid f 


A 70 -1- fl fl/l 
4. / z ± U.U4 


a ox. -i- n n/i 

4.»0 ± U.U4 


c OK -I- n OK 

o.yo ± u.uo 


n oc; 4- n no 
u.yo ± U.Uz 


n oo 4- n m 
u.y^ ± U.U1 




vise fluid 2 


4.oU ± U.U4 


/i no -J- n n/i 
4.yy ± u.U4 


^ OQ -J- O OK 

o.yo ± u.uo 


n ok 4- n no 

u.yo ± u.uz 


U.yo ± U.U1 




vise fluid 3 


A co -L n r\A 

4.0Z ± U.U4 


a ok -L n n^ 

4.yo ± u.uo 


p. on J- o ok 
D.UU ± U.UD 


n o7 4- n no 

u.y/ ± u.uz 


n oo i n m 

u.y^ ± u.ui 




vise fluid 4 


4.oy ± U.U4 


k no 4- n n^ 
O.Uz ± U.UO 


p. 1 K _l_ O OK 

D.1D ± U.UD 


n 07 4- n no 
u.y/ ± U.Uz 


O OQ 4- o 01 

u.yo ± u.ui 


0.39 


ideal fluid 


4.68 ± 0.06 


4.73 ± 0.07 


5.06 ± 0.08 


0.99 ± 0.03 


0.93 ± 0.02 




vise fluid 1 


4.63 ± 0.06 


4.77 ± 0.07 


4.91 ± 0.08 


0.97 ± 0.03 


0.92 ± 0.02 




vise fluid 2 


4.66 ± 0.06 


4.76 ± 0.06 


4.85 ± 0.07 


0.98 ± 0.03 


0.91 ± 0.02 




vise fluid 3 


4.67 ± 0.06 


4.79 ± 0.07 


4.95 ± 0.08 


0.97 ± 0.03 


0.91 ± 0.02 




vise fluid 4 


4.77 ± 0.07 


4.81 ± 0.07 


5.09 ± 0.08 


0.99 ± 0.03 


0.90 ± 0.02 


0.55 


ideal fluid 


4.46 ± 0.08 


4.51 ± 0.08 


4.38 ± 0.08 


0.99 ± 0.03 


0.95 ± 0.02 



vise fluid 1 4.45 ± 0.08 4.48 ± 0.08 4.28 ± 0.08 0.99 ± 0.04 0.95 ± 0.02 

vise fluid 2 4.49 ± 0.08 4.52 ± 0.08 4.17 ± 0.08 0.99 ± 0.04 0.94 ± 0.02 

vise fluid 3 4.49 ± 0.08 4.63 ± 0.09 4.40 ± 0.09 0.97 ± 0.04 0.90 ± 0.03 

vise fluid 4 4.52 ± 0.08 4.60 ± 0.09 4.48 ± 0.09 0.98 ± 0.04 0.92 ± 0.03 



0.70 


ideal fluid 


4.25 ± 0.09 


4.51 ± 0.10 


3.88 ± 0.09 


0.94 ± 0.04 


0.98 ± 0.03 




vise fluid 1 


4.23 ± 0.09 


4.56 ± 0.10 


3.92 ± 0.09 


0.93 ± 0.04 


0.99 ± 0.03 




vise fluid 2 


4.23 ± 0.09 


4.59 ± 0.10 


3.78 ± 0.09 


0.92 ± 0.04 


0.95 ± 0.03 




vise fluid 3 


4.26 ± 0.10 


4.56 ± 0.11 


3.91 ± 0.10 


0.93 ± 0.04 


0.95 ± 0.03 




vise fluid 4 


4.30 ± 0.10 


4.67 ± 0.11 


3.89 ± 0.10 


0.92 ± 0.04 


0.94 ± 0.03 


0.84 


ideal fluid 


4.23 ± 0.12 


4.42 ± 0.13 


3.37 ± 0.10 


0.96 ± 0.05 


0.97 ± 0.04 




vise fluid 1 


4.24 ± 0.12 


4.46 ± 0.12 


3.43 ± 0.10 


0.95 ± 0.05 


0.97 ± 0.04 




vise fluid 2 


4.32 ± 0.12 


4.40 ± 0.12 


3.59 ± 0.10 


0.98 ± 0.05 


1.01 ± 0.04 




vise fluid 3 


4.44 ± 0.13 


4.48 ± 0.13 


3.60 ± 0.11 


0.99 ± 0.06 


1.02 ± 0.04 




vise fluid 4 


4.77 ± 0.14 


4.85 ± 0.15 


3.67 ± 0.12 


0.98 ± 0.06 


1.05 ± 0.05 
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FIG. 6: (Color online) The two-pion HBT results for the viscous (C s = 2) and ideal hydrodynamic 
sources with eq = 20 GeV/fm 3 and eo = 40 GeV/fm 3 . The freeze-out temperature is 0.8T C . The 
experimental data are from PHENIX [47[ and STAR [48j] Collaborations. 



onance gas at low temperatures with lattice QCD at high temperatures 44j. We consider 
both the shear and bulk viscosities in the hydrodynamic evolution of the dissipative fluids. 
By performing the two-pion HBT interferometry to the viscous and ideal hydrodynamic 
sources, we investigate the effects of the viscosities on the two-pion HBT interferometry 
results. 

For a dissipative fluid there is not only shear viscosity but also bulk viscosity. For the 
sQGP produced in relativistic heavy ion collisions [lj], the ratio of the shear viscosity to 
the entropy density satisfies rj/s >l/47r [38l-l40| . and the bulk viscosity ( is proportional to 
r)(Vg — 1/3) 18|, l4ll443j]. Our calculations indicate that the viscosity effect on the two-pion 
HBT results is small if only the shear viscosity is considered. The bulk viscosity leads to 
a larger transverse freeze-out configuration of the pion sources, and therefore increases the 
transverse HBT radii. This effect is larger if the source has a lower freeze-out temperature. 



In the present paper we took the relaxation time tv 



6)7 
sT 



13 



16 



and r n = r n 
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, and used an uniform formula of the shear viscosities for both the QGP and hadronic 
phases for simplicity. We find that there are not observable influences on the viscosity 
effects on HBT results for the r n and rn values within reasonable ranges. In Ref. [19 1 
different viscosities for the QGP and hadronic phases are taken into account in the analyses of 
transverse momentum spectra, elliptic flow, and HBT interferometry. Further investigations 
of the effects of shear and bulk viscosities on the HBT interferometry for more realistic 
viscous pion-emitting sources in relativistic heavy ion collisions are of great interest in the 
future. 



In our model the Bjorken longitudinal scaling hypothesis 26] is adopted, which is an ap- 
proximation for the heavy ion collisions at very high energies. Our results indicate that the 
longitudinal HBT radius for the Bjorken cylinder hydrodynamic source agrees well with the 
RHIC HBT results. However the transverse HBT radii as functions of the pair transverse 
momentum for our Bjorken cylinder sources are more flat than those of the experimental 
data. In order to explain the RHIC HBT results, various models without viscosity were pro- 



posed 32|, l49H54j . Recently, the hydrodynamic source model including the initial collective 
flow, stiffer EOS, and viscosity was proposed |55j to explain the RHIC HBT results. The 
comprehensive examinations of the effects in these models on HBT interferometry and other 
observables (e.g. particle momentum distributions, elliptic flow, etc.) require considerable 
efforts. 
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Appendix A: Derivation for relaxation equations 



Here we give the derivations of Eqs. (1241) - (I26p from Eq. (I22p in detail. For the Bjorken 
cylinder and at z = 0, we have 
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and 



4 = gT^ 



(A3) 



(l 
-1 

\0 

/ ^ 2 v 2 7 pp 
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1/7* J 




o -r 
o o 



-"f 2 V p T pf> 



-7 2 V p 7 pp 
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-7 2 V p 7 pp 
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T 2 T VV J 



(A4) 



1. The first term of Eq. ([22]) 



The first term of Eq. fl22l) can be expanded as 

- t^u^Dtt^ + 7^u a u v u p Dw a p. 



(A5) 



With the relations Dg K p = Dg Kl3 = and the orthogonality 



0, we have 



A^Dv = Al(^) - (^) ^] = -<TV(W) <n a p 
g^u^Dir^ = g^ U P[(Du a n^) - (Du a ) n a0 ] = -g^u^Du") tt q/3 
u / Vu £ Vr>7r a/3 = uV[u a (Du P 7r a0 ) - (Du 13 ) u a 7r al3 } = . 
So, the first term of Eq. f|22|) can be expressed as 



where 



= <TV(Di/) vr Q/3 + g^(Du a ) vr Q/3 . 



(A6) 

(A7) 
(A8) 
(A9) 

(A10) 
(All) 



For the azimuth-symmetric Bjorken cylinder, = and u v = at z = 0, so 1^ = I m = 0. 
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Therefore, for p = v = p we have 



A^A^Dn^ = g pl3 Dn p + I pp = D-f 2 r pp + I pp 
d_ 
'dp 



27(i57) T PP + I PP , 



(A12) 



where 



Ipp = g pa u p {Du^)Ti a p + g pp u p {Du a )ix a p = g pa u p {Du^)Tx a p + g pa u p {Du p )-K Pa 
= 2g pa u p (DuP) 7i aP = 2g pp u p {DuP)<K p p = -2^[( J D 7 )vr /9T + (D<yv 
= -2jv p [- (Dj)j 2 v p t pp + (Djv p )j 2 t pp ] = -2j\t pp (Dv p ). (A13) 



With the relations, 



7 



<9( 7 f p ) _ ®1 d^yvp) _ dj 



y/l — v 2 p ' Up dr dr' Vp dp dp 



(A14) 



we have 



pp = -2^T pp [v p (D 1Vp ) - v 2 (D^)] = -2 7 V P 



-2 7 V 



dr 

-2jt pp Dj, 



and 



d 



-pp 



dp 



for p — v — p. 

For p = v = cf) and p = v = 77, we have 



A^A"? DlTaf, = —Dt^ 



p 2 \dr 



d_ 

'dp 



and 



1 f) 



(A15) 



(A16) 



(A17) 



(A18) 



2. The third term of Eq. rt22l 



For the third term of Eq. (1221 , we have 

V <p u u> = V {p u v) - -A pv d x u } 

3 



W + V y u p 



-Are. 
3 



(A19) 



20 



With the relations 



we have 



and 



' dr dr ' dp v p dp 

dj d(^v p )- 



V p u p 
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-dr dp 



-Ye. 



(A20) 



(A21) 



where 



Similarly, we have 



V <p u 
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dp 
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■■>..<•>- 1 f-^p | 1 e 
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T V T 
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p w 5 



V <T ?i^ = V <p u T> = 7 2 v p ( - + _e) = 7 2 w p a Tp 



(A22) 



(A23) 

(A24) 
(A25) 
(A26) 
(A27) 

(A28) 

(A29) 
(A30) 
(A31) 



3. The last term of Eq. ([22 



The last term of Eq (l22|) can be written as 



//■Vv/ A ( ) V 7 " 



7/ i \n 1 / \n 1 



-r]T 7 



^p T/ 



(A32) 
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Now, one can obtain Eqs. (|24|) - (I26I) from the above derivations. 
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